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This paper presents a comparison of the theoretical and experimental current-voltage (/-V) 
characteristics of a self-assembled monolayer of ar,a-xylyl dithiol molecules on a gold substrate 
measured with a scanning tunneling microscope probe. Good quantitative agreement is obtained with 
the tip-molecule distance as the only "fitting parameter** Several other thiol<ouplcd molecules 
that we have studied also show similar agreement The conceptual picture presented in this paper 
could be usefiil for the interpretation of I-V measurements on molecular monolayers in general. 
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Scanning tunneling microscope (STM) images of self- 
assembled monolayers of different molecules on different 
substrates (typically gold, [1]) have been reported by many 
groups. Such images arc obtained by recording the z co- 
ordinate of the tip as it scans the surface at a fixed current 
[2]. However, there are relatively few reports of current- 
voltage (/- V) characteristics through such molecular mono- 
layers [3-6]. Quantitative comparisons of the theoretical 
and experimental /-K characteristics are limited to the low 
bias regime (for Ceo molecules [7]), or to structures with 
intervening metallic clusters whose I-V is limited by the 
Coulomb blockade due to the cluster [8,9]. This paper 
presents what is probably the first report describing a quan- 
titative comparison of the theoretical and experimental I-V 
characteristics of individual molecules. Good agreement 
is obtained with the tip-molecule distance as the only "'fit- 
ting parameter," providing support for the basic picture for 
molecular conduction reported last year by a number of 
groups [ 1 0- 1 3]. However, this agreement is obtained only 
if the potential drop between the molecule and the substrate 
is taken into account. The conceptual picture presented in 
this paper could be useful for the interpretation of /- V mea- 
surements on molecular monolayers in general. 

Molecular monolayers generally attach strongly to the 
gold surface through thiol ( — S — ) end groups [14.15] giv- 
ing a strong chemical bond with good orbital overlap. The 
other contact, typically an STM probe (see Fig. 1), is usu- 
ally coupled weakly to the molecule. It is commonly as- 
sumed that the applied potential eV is dropped entirely 
between the tip and the molecule, while the molecular 
potential remains fixed with respect to the substrate [see 
Fig. 2(b), with V^oi 0], For a positive substrate voltage 
(as shown) the molecule conducts str ngly as the electro- 
chemical potential in the tip. approaches the lowest un- 



occupied molecular orbital (LUMO); for negative substrate 
voltage the molecule conducts strongly as /12 approaches 
the highest occupied molecular orbital (HOMO). Based on 
this picture we should expect the /-K characteristics to be 
strongly asymmetric for all thiol-coupled molecules since 
the HOMO is a sulfur-based level that couples strongly to 
the gold while the LUMO is a ring-based level that couples 
only weakly. The current should thus be much smaller for 
positive substrate voltage (conduction through the LUMO) 
than for negative substrate voltage (conduction through the 
HOMO). However, all the molecules we have studied so 
far show symmetric I-V characteristics (unless the tip is 
moved very far from the molecule). We see no simple ex- 
planation for this observation as long as we assume that 
the molecular potential remains fixed with respect to the 
substrate. But. as we will show in this paper, the symme- 
try of the I-V characteristics is explained quite easily once 
we take the shift in the molecular energy levels under bias 
into account. 




FIG. 1. Self-assembled monolayer of a, a -xylyl dithiol on a 
gold substrate (only one molecule of the monolayer is shown 
for clarity). Also shown schematically is an STM probe used 
to measure the current-voltage characierisiics. 
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IVkat is the potential profile? — We assume that the 
current flows through an individual molecule without 
any significant lateral conduction through the molecular 
monolayer. The assumption that the molecular p tential 



is fixed with respect to the substrate is usually justified 
by arguing that the substrate-molecule resistance R\ is 
much smaller than the tip-molecule resistance /?2. Since 
the current / is the same across either junction, one can 
write [see Fig. 2(a)] 



Mmol - Ml M2 - Mmol Mmol " Ml 



/?2 



as /?2/^i °°- However, what this argument proves 
is that the electrochemical potential in the molecule 
(Mmoi) >s the same as that of the substrate (^ti), as 
shown in Fig. 2(a) [16]. It says nothing about the 
electrostatic potential (p( ) which can vary spatially in 
inhomogeneous conductors even in equilibrium. Such 
equilibrium variations in <p( ) are already included in the 
calculation of the molecular energy levels. What we need 
to know is the change 8<p{ ) in the electrostatic potential 
under bias which enters the molecular Hamiltonian and 
affects the molecular energy levels (the electrochemical 
potential enters only indirectly through charging effects). 

To see what 5<p{ ) looks like, we first note that inside 
the substrate or the tip it must be zero, since these regions 
are charge neutral. If we neglect any charge buildup inside 
the molecule, S<p{ ) inside the molecule can be obtained 




TransmtMlon TIB) 



FIG. 2. (a) If we assume that the electrostatic potential 
varies linearly from the substrate to the tip (C), the average 
electrostatic potential in the molecule (Vmoi) will be about half 
the applied voltage V. The electrochemical potential (^4) in the 
molecule is approximately the same as that of the substrate, (b) 
Transmission ftmction T{E) through the molecule fiom the tip 
to the substrate with (solid) and without (dashed) scattering. 
Taking the molecular potential as our reference, ^\ moves 
down by eV^i while fxz moves up by ^IV-Vntoi) where V^^i is 
the average electrostatic potential in the molecule. 



by solving the Laplace equation V^(5^) 0. If we view 
the STM tip and the gold substrate as the two plates of 
a parallel plate capacitor then the electrostatic potential 
will vary linearly, as shown in profile C of Fig. 2(a), 
and the average electrostatic potential in the molecule 
Vnioi can be written as Vmoi/ V Cmoi/^. where Znoi and 
L are the distances from the substrate to the center of 
the molecule and the tip, respectively. To the lowest 
order of approximation, we can assume that the molecular 
energy levels simply float up by an amount eVmoi- More 
conveniently, we can take the molecular energy levels as a 
fixed reference and let the substrate float down by eVmoi* 
as shown in Fig. 2(b): 



/ii Ef - T)eV 



and 



M2 + v)eV (I) 

(note that e -1.6 X 10"*^ C), where the factor rj de- 
scribes how the electrostatic potential difference V is di- 
vided between the two junctions: 17 = Vmoi/V Zmo\/L. 
Since the tip has to be within a few angstroms firom the 
end of the molecule in order for the current to be mea- 
surable, Vnjo! is approximately half the applied voltage V 
in = 0.5). 

Effect of 7) on the I-V characteristics.— Fig\m 2(b) 
shows the transmission ftmction through the molecule 
fi'om the tip to the substrate which is calculated using 
an extended Huckel model for the molecule and taking 
the metal-molecule bonding into account. It shows peaks 
at energies corresponding to the molecular energy levels 
which are broadened further if we include scattering 
processes in the molecule. The strong asymmetry in the 
transmission for the levels above and below the gap is 
evident. Depending on the value of 77 (that is, V„,oi), 
we can obtain very different predictions for the I-V 
characteristics. If 77 0, /xi is fixed with respect to the 
molecule. The molecule then conducts strongly when ^2 
coincides with the energy of a molecular orbital, that is, 
when 



or 



eV > El - Ef (positive substrate voltage) 



—eV > Ef - Eh (negative substrate voltage). 



The conductance spectrum shows a gap of (E^ — Eh) 
equal to the molecular gap. This is the usual case in 
STM measurements. But consider what happens if 17 is 
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approximately 0.5 and neither /xi nor fi2 remains fixed 
with respect to the molecule. The molecule then conducts 
strongly when either /xj or coincides with a molecular 
orbital energy so that the threshold for conduction is 
given by 

. ( Ef - Eh £ljzIl\ 

eV > mml , I 

\ 77 1 - 77 / 

(positive substrate voltage) (2a) 

or 

. ( Ef - Eh Et-Ef \ 

V 1 - 77 T) J 

(negative substrate voltage) , (2b) 

where min(jcj') denotes the smaller of the quantities x 
and y. 

It is easy to see from Eq. (2) that with 77 = 0.5 the 
molecule could conduct through the HOMO for both bias 
polarities leading to symmetric I-V characteristics, if the 
equilibrium Fermi energy £/ is closer to the HOMO than 
to the LUMO. This is what we find in all the thiol- 
coupled molecules when we locate the equilibrium Fermi 
energy £/ using the relation 

N 2(for spin) X Y — tan" * — , (3) 

where F, is the broadening of the energy level located at 
e, and N is the number of electrons in the molecule [17]. 
The procedure used to calculate f/, is described below. 
For a, a -xylyl dithiol [shown in Fig. 1(a)], we find 

Ef - Eh 1.8eV, El ' Ef 2.2eV. (4) 

From Eqs. (2) and (4) we expect the threshold bias 
for conduction to be approximately 3.6 V for either bias 
polarity if we assimie 77 0.5. The calculated /-F based 
on this assumption fits the experimental data quite well, 
as shown in Fig. 3(a). A calculation with 77 0, on the 
other hand, shows hardly any current for positive substrate 
bias and is clearly in disagreement with experiment. 

Theory. — Let us briefly describe the procedure used to 
calculate the I-V shown in Fig. 3. We use the standard 
expression for the current based on the scattering theory of 
transport [18]: 

/ Y f ^ dET(E)[f(E - /i,) - /(£ - PL2)]. 

(5) 

Here T{E) is the transmission function from the tip to the 
substrate and /(£) is the Fermi function. The electro- 
chemical potentials /xj and ^2 are obtained from Eq. (1), 
with the equilibrium Fermi energy £/ located from Eq. (3) 
as described earlier. We use the extended Hiickel method 
[19] to obtain the Hamiltonian (//) for the isolated mole- 
cule and include the effect of the contacts through self- 
energy functions (S) that are calculated taking the proper 
metal-molecule coupling into account [20]. The eigen- 
energies of (// + 2) are complex: £, e, + /T,. The 




Substrate Voltage (V) 
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FIG. 3. (a) Current- voltage characteristic for the molecule 
shown in Fig. 1 plotted on a linear scale: 77 0.5 (solid), 77 

0 ( )» and measured (+). (b) Current- voltage characteristic 

for the molecule shown in Fig. 1 plotted on a logarithmic scale: 
77 0.5 (solid), and measured (+). Horizontal dash-dotted 
line represents the noise limit of our preamplifier. 

real part (e, ) is the energy of a level while the imaginary 
part (F,) represents its broadening due to the coupling to 
the metal and describes how effectively the level can empty 
into the metallic reservoirs. The transmission function 
r(£) is calculated ft-om H and S as described in [12]. It 
shows peaks at energies corresponding to the molecular 
energy levels [see Fig. 2(b)]. If we include scattering 
processes in the molecule through an additional self-energy 
function [21], the peaks are broadened further, as we would 
expect. This has no significant effect on the /-F character- 
istics shown in Fig. 3, apart from rounding out the cusps. 

Results — Figure 3(a) shows the experimental I-V 
characteristics measured in ultrahigh vacuum (-^3 X 
10"*° Torr) with a set current of 0.25 nA at 5 V. The 
experimental procedure is essentially the same as that 
described in Ref. [8], except that there is no intervening 
cluster between the tip and the molecule. There is good 
agreement between the experiment and the theoretical 
calculation from Eq. (5) with 17 0.5. By contrast, with 
77 0 there is hardly any current for positive substrate 
bias due to the weak coupling to the LUMO, in clear 
disagreement with experiment. Figure 3(b) shows the 
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/-K characteristics on a logarithmic scale showing that the 
agreement (with 7; 0.5) is good over several orders of 
magnitude of the current. Note that there is no fittmg 
parameter in this calculation apart from the distance from 
the tip atom to the sulfur atom, which is adjusted to give 
the same magnitude of the current as the experiment [22]. 

Discussion and conclusions.— Ix is evident from Fig. 3 
that there is clear disagreement with experiment if we as- 
sume that the electrostatic potential follows curve A m 
Fig 2(a) (the same as the electrochemical potential). We 
see no simple way to explain the observed characteristics 
based on profile A. But excellent agreement is obtamed if 
we assume that the electrostatic potential follows profile B 
with Vmoi 0-5 V (that is, v 0-5)- The actual electro- 
static potential should be calculated from the Poisson equa- 
tion V • (eVS(p) Bp, where the change in the charge 
density (5p) has to be determined self-consistently [23]. 
An exchange-correlation potential corresponding to 5p 
should also be included in the calculation of the molecu- 
lar energy levels. We find 6 p to be negligible for most 
of the bias range considered here. With 5p = 0, the po- 
tential profile would be linear (profile Q. However, the 
agreement with experiment deteriorates significantly if we 
use profile C instead of profile 5, leading us to believe 
that the latter is closer to the actual profile. This could 
be justified by arguing that the molecule has a higher di- 
electric constant e compared to the gaps at the ends. The 
theoretical model presented here should be viewed only as 
a first step towards a more complete theory that mcludes 
other details such as the correct self-consistent potential 
profile or the structure in the density of states in the con- 
tacts. However, we find that the simple model presented 
here leads to good agreement with experiment not only for 
a a -xylyl dithiol (described in this paper), but also for 
a number of other thiol-couplcd molecules that we have 
studied. In some cases, an adjustment of the eqmlibnim 
Fermi energy £/ leads to better fits with experiment; this 
can be justified by the fact that £/ is very sensitive to small 
amounts of broadening and charge transfer. Also, the I-V 
characteristics become asymmetric (like the ones reported 
in Ref. [4]) when the tip is moved very far from the mole- 
cule, but the data can still be described by our model using 
a suitably modified t/. 
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